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Low cycle fatigue has been a subject of extensive study [1,2]. Fatigue 
damage consists of crack nucleation and crack propagation. A small crack can be 
nucleated under a cyclic load by repeated cyclic slip, creep void nucleation and 
void coalescence at elevated temperatures, oxidation or other environmental at- 
tacks on grain boundary, fracture of oxides, intermetallic compounds or other 
inclusions, and debonding at the particle-matrix interface. Once a micro-crack 
is nucleated, and the crack is large enough, it will propagate to failure under 
a continued cyclic load. Therefore, the total fatigue life, N^, consists of two 
parts: N Q and Np for crack nucleation and propagation albeit the dividing line 

between nucleation and propagation is not definitive. 

It is generally agreed, at the higher applied load range, that the relative 
portion of the nucleation life to the total life, N Q /Nf» is lower. Therefore, in 
the low cycle fatigue range, fatigue life is consumed mainly for crack propagation, 
especially where environmental attack is severe. 

It is necessary to extend fatigue crack growth analysis into the region of 
general yielding in order to study low cycle fatigue life. This paper analyzes 
the crack tip field in general yielding, relates the crack tip field to far field 
parameters, analyzes Solomon's fatigue crack growth rate measurement [3] in terms 
of characteristic crack tip field, i.e. in terms of J, develops a general expression 
for fatigue crack growth, derives a correlation for low cycle fatigue life and com- 
pares the derived correlation with experimental data. 
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I. CRACK TIP FIELD IN GENERAL YIELDING AND FAR FIELD PARAMETERS 

Hutchinson [4] and Rice and Rosengren [5] studied the characteristic crack 
tip field in an elastic-plastic solid. According to their analyses, the crack 
tip field can be characterized by a single parameter, K o , or K^, or J, as defined 
by Hutchinson, and Rice and Rosengren. The crack tip paramters, J, K and K are 

O £ 

related to the far field parameters, and the detailed relation depends on specimen 
geometry as well as on the stress-strain relation. J is defined as a contour in- 
tegral, and it can be determined experimentally for elastic-plastic solids by the 
deformation work difference AWp as shown in Fig. 1.1. 

Under a cyclic loading, the stress-strain relation in the plastic region often 
exhibits a power relation. Three different power hardening relations have 
been used to simulate the real behaviors of materials: 

1. Pure power hardening, pure PH 



where, o Q , e Q , are reference stress and strain, and a and N are constants. 

2. Ramberg-Osgood power hardening - Ramberg-Osgood PH 

£_=2_ +a (E_) N (1.2) 

e o °o °o 

with o Q = Ee 0 . 

3. Piecewise power hardening - piecewise PH 

a = Ee o <_ Oy (1 - 3a) 

— = (— ) n o > o Y (1.3b) 

Oy 'Ey' Y 

The stress-strain relations of piecewise power hardening materials can also 


be simulated by 

■ . - + C '■ °V 

£y Oy Oy 

>- 

D 

> 

A 

D 

(1.3c) 


o = Ee 

O < Oy 

(1.3d) 
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The second term of the right hand side of Eqn. (1.3c) is the plastic strain. 

It is applicable only when o > Oy. 

In this section, the earlier development of the characteristic crack tip 
field in large scale yielding and general yielding will be reviewed. The de- 
tailed characteristic crack tip field was derived by Hutchinson and Rice and 
Rosengren for pure PH materials. The relations between J and far field para- 
meters for pure PH materials have been studied by Goldman and Hutchinson [6] 
and Shih and Hutchinson[7] for plane strain and plane stress cracked panels in 
tension. Amazigo [8] analyzed the pure PH solid in anti-plane shear. Rice [9] 
related J to the applied deformation work density for a piecewise PH semi- 
infinite-cracked infinite-strip under antiplane shear. Shih [10] obtained the 
antiplane shear results for both piecewise PH and Ramberg-Osgood PH materials 
for a finite specimen. Shih and Hutchinson [7] have also studied the plane stres§ 
center-cracked panel in tension for Ramberg-Osgood PH solid. 

In this section, the relation between J and the far field parameters, the 
applied stress, o b , the applied strain, and the deformation work density, Wp 
are studied for piecewise PH materials by FEM. The pure PH relation neglects the 
elastic part of the stress-strain relation entirely. For Ramberg-Osgood PH re- 
lation, J is related to two terms of far field parameters [7,10]. For piecewise 
PH relation, J is related to only one term of far field parameter, which is 
simpler in analyzing low cycle fatigue behavior. 
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1.1 Characteristic Crack Tip Field- Elastic-Plastic Solids 

Goldman and Hutchinson [6] formulated the fully plastic crack problems to 
deduce the crack tip field characteristics. With the assumptions of (1) small 
strain deformation theory, (2) incompressible solids, (3) proportional loading 
and (4) pure power hardening stress-strain relation, Ilyushin concluded that if 
the boundary traction increases linearly in proportion to a single load para- 
meter t, the stresses at each point in the solid also vary linearly with t 
[11,6,12]. With the pure power hardening relation. 


or 


where 



N 


11 = T °<^> 


- N -’ ill 


1J 


- o 


1J 


3 °pp 6 ij and °e " 2 S ij S ij 


(1.4) 


(1.5) 


the stresses, strains, and displacements of the characteristic crack tip field have 
the form [13]: 
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u. = 
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ae K 
0 E 


r 1+N ^(e.N) 


where r and e are the polar coordinates with the origin at the crack tip and the 
crack line lies along the line e = n. K and K are stress and strain intensity 

m O G 

factors, and they are related by (K ) = K . a..(e,N), o.(e,N), e..(e,N), and 

'Y, O C 1 J 6 1J 

u^e.N) are the dimensionless functions of e and N. They define the e-distributions 
of their corresponding stress, strain or displacement components. a Q and e Q are 
reference stress and strain, a and N are the dimensionless constants. 


4 


The strain energy of a non-linear elastic Ilyushin solid is defined as 


W = 


f °ij 


j _ /_N_w_e\l+N 

°kl dc kl " a °o e o'l+N o Q 


(1.7) 


Rice [9] defined J-integral as 


J = 


(Wn, - o j:j Hj u (> 1 ) ds 


0 . 8 ) 


and J is related to K and K 

0 E 


J = 


o a e K K I = 
o o a e n 


(Y T 

a o e K I 

o o a n 


1+N 

a o e (K ) ^ I 
o o' e n 


(1.9) 


where I is a numerical constant dependent on N, and the values are given in 
Ref. [13]. 

Equations (1.6) define the crack tip field, and they are asymptotic solutions 
valid only when r approaches the crack tip. These equations give the characteristic 
crack tip field for non-linear elastic Ilyushin solids. It has been accepted that 
these equations also approximate the crack tip field in power hardening elastic- 
plastic solids; and Equations (1.6) have been widely used in elastic-plastic frac- 
ture mechanics. Indeed, if there exists a unique crack tip field invariant to in- 
plane geometric variations and the types of loading, these equations together with 
the relations for J, K o and K £ will be able to form a basis for a generalized 
elastic-plastic fracture mechanics theory. However, McMeeking and Parks [14] and 
Shi h and German [15] have shown that the plane strain crack tip field is dependent 
on sample geometry and the types of loading. Therefore, the characteristics of the 
crack tip field in general yielding have to be examined carefully before the use of 
a mechanical parameter such as J, K o> K^, and rp to characterize a specific crack 
tip field. 

To begin with, a reference characteristic crack tip field has to be chosen. 

The elastic-plastic crack tip field in the condition of small scale yielding seems 
to be the best choice as the reference field. Hu [16] and Hu and Liu [17] have 
suggested the use of direct correspondence between the crack tip field in a large 
sample in small scale yielding and the crack tip field in a small sample in general 
yielding as the basis of elastic-plastic fracture mechanics. 
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1.2 J and Far Field Parameters 

For Ilyushin solids, the crack tip stresses and strains are linearly propor- 
tional to the boundary traction, o . Therefore, the crack tip field can be ex- 
pressed in the form [6] 




1 a 

a ’ w 


.N) 


( 1 . 10 ) 


and 


1 j _ ^ /X y a 

7 *" " o(— ) E..(- , t . »N) 
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J = a a e a 

o o o 


(-V •"> 


( 1 . 11 ) 


where a and w are crack length and panel width, and the dimensionless functions 
topped by 'V are independent of the boundary stress o o . 

Plane stress and plane strain crack tip fields in pure PH materials have been 
calculated by Goldman and Hutchinson [6] and Shih and Hutchinson [7] for center 
cracked panels in tension. J is related to the applied stresses, and the relations 
are given in the following form: 


J 

ao e a 
o o 



for plane stress 


( 1 . 12 ) 
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for plane strain 


(1.13) 


The values of g and f for pure PH materials are tabulated over the ranges 
g <J<| and 1 < N < 20 [6,7]. 

For both Ramberg-Osgood PH and piecewise PH materials, the stress-strain re- 
lations do not fulfill the requirements of an Ilyushin solid. Therefore Eqn. (1.11) 
has to be modified. 

Rice [9] has shown that the J-integral of a semi-infinite crack in an infinite 
slab with clamped boundaries at a constant shearing displacement is 

J = 2hW 

oo 

where 2h is the height of the slab. In an elastic solid, is the constant 
strain energy at x = + «*> as shown in Fig. 1.2 . For piecewise PH materials. 


6 


(1.14) 
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J is given by [9] 
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(1.15) 


Shi h and Hutchinson [7] have obtained the numerical results for plane-stress 
center-cracked panel in tension for Ramberg-Osgood PH materials. This section 
deals with the numerical calculations of plane stress and plane strain single- 
edge cracked panels in tension using piecewise PH relation following the procedure 
of Shih and Hutchinson [7]. 

For piecewise PH materials in tension 


o 




a £ 0y 



a > Oy 


(1.16) 


where Oy and Cy are yield stress and strain respectively. 

Based on Eqns. (1.12), (1.13) and (1.15), the expression for J-integral can 
be generalized to the following form: 
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(1.17a) 
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(1.17b) 
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where B = 1 for plane stress case, and B = (jp) n for plane strain case. 

9 is a function of (^) and n. The values of 9 for plane stress and plane strain 
cases are different, and they are tabulated in References [6,7]. As indicated 
by Goldman and Hutchinson [6], for plane strain case, when x 2 approaches «, 
we have o n - 0, a 22 - o^, a 33 - j a*. -> ^ o m ,e 33 = 0 , c 22 = - . 

From Eqn. (1.5), we also have the following for the pure PH materials in the plane 
strain case: 


e °° _ /3 / /3 g <»\n 

17 ' 2 l 2 y 


0-18) 


and 


a e 1+ - 

°° °° _ / /J «” , n 

0yEy 2 0y 


(1.19) 


In Eqn. (1.17b), the first term is the elastic contribution, and the second 
terms are the plastic part. 

The expressions for the stress intensity factors in the linear elastic 
fracture mechanics are modified by using a g yy in order to take account of the 
plastic effects [18]. To account for this effect, following the practice by 
Shi h and Hutchinson for the Ramberg-Osgood PH materials[7] , Eqns. (1.17a and b) 
are modified to the following forms for piecewise PH materials. 
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where 


a eff = a + r p 
a eff = ^eff^a = 0 , 


for 0 < o v 
00 — y 


for 0 > 0w 

* y 


( 1 . 21 ) 


K t 2 
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for plane stress 


for plane strain 


( 1 . 22 ) 


eff /W-a \ 

1 ~^ w - a eff 


(1.23) 


For the study of fatigue lives of smooth specimens, the relevant fatigue 
crack is often very small relative to the other dimensions of a specimen. There- 
fore, the case - -*■ 0 is of particular importance. In this case: 
w 


, . a eff, w-a , a eff r , , 1 ,1-n w°»> 2 i 
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g(~^, n=l) = g(0,l) 
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I 00 


Equations (1.20a and 1.20b) become 
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0„£0 y (1.24a) 
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where g(0,n) is given by 

g{0,n) = [3.85/1 (1-n) + Tin] 

for a center cracked panel in the plane stress case [7]. A similar formula 
can be written for the plane strain case. 

Equation (1.24) can also be written as 
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Wp is far field deformation work density as illustrated in Fig. 1.3, and it is 


given by Wp 


o de . The size of a critical fatigue crack at failure is 

GO 00 

Jo 


usually rather small. Therefore, for a very small crack, o o and e o approach o 
and e of a smooth specimen. 

The deformation work done to a non-linear elastic solid is stored as poten- 
tial energy , i.e. strain energy. But for an elastic-plastic solid, a large 
part of the deformation work is dissipated in the form of heat; only a small part 
is stored as potential energy. 

Equations (1.24) and (1.25) can be written in the simpler form: 


j = g(0,l)(o„cja eff 

= CnJ + n 2 o y e Y ]a , 

and 

J = 2g(0,l)Wp a eff 

II II 

= [n 1 Wp + n 2 a y c y ]a , 

I I ll ll 

where n-| t n 2 » n-j and n 2 are functions of n. For plane strain case, the re- 
lation between J and the far field parameters can be expressed in the same form. 
Equations (1.27a and 1.27b) will be used as the basic forms of the expressions 
of J for a small crack in a wide plate of a piecewise PH material. These two 
relations will be checked with the finite element calculations and the values 

I I II II 

of rj -j , n 2 » n-j , and n 2 will be determined. 

1.3 The Finite Element Method Analysis 

In order to- assess the validity of the above analysis, an elastic-plastic 
finite element study was undertaken. The ABAQUS FEM program Was used. The 
finite element idealization of a single edge cracked panel in tension is shown 
in Fig. 1.4. The crack length to width ration, a/w = 0.1. The material obeyed 
the Von Mises criterion and the piecewise power hardening relation. 
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o = Ee ; 


a — a y 

f- = (|-)" S o > o ¥ 0.28) 

The exponent of the stress-strain relation, n, was taken 0.2 and 0.3. 

4 

E = 7.0 x 10 MPa and a y = 350 MPa were assumed. 

The plane eight noded quadratic isoparametic element was used. The crack 
tip triangular quadratic isoparametric elements were formed by collapsing one 
side of the quadratic element and placing the mid-side node closer to the crack 
tip at the quarter point [193 - Crack-tip blunting is modelled by separating the 
condensed crack tip nodes as shown in Fig. 1.5. Displacement control was used. 

Both plane stress and plane strain analysis were made. 

J-integral are evaluated from the field values obtained by the finite ele- 
ment analysis. Both rectangular path and circular path are chosen. The paths 
are shown as dashed lines in Fig. 1.4. The detail procedures of the calculation 
were reported by Zhuang [20]. 

Path dependency are less than 3% for all the plane stress cases and less than 
7% for all the plane strain cases. 

Figures 1.6 to 1.13 show the plots of o,„/o v , e w „/e w , e P /e v , and o /o v along 

yy y yy y e y eY 

the crack line versus r/(d/o y). J/<jy is a measure of crack opening displacement. 

In Figs. 1.6 to 1.9, the different symbols represent the data at different boundary 
displacement levels. The solid symbols indicate that the applied nominal strain 
is beyond the yield strain, ey. The open symbols indicate that the nominal strains 
are below ey. The nominal strain is the quotient of loading displacement divided 
by height of the specimen. In each dimensionless plot, all the data near crack 
tip for small scale yielding, large scale yielding and general yielding fall on 
the same curve. The fact that, at various load levels, the data of e yy , a g , 
and e p , in the crack tip region fall on the same curves suggest that there exists 
a characteristic crack tip field and that a single parameter J, or can 
characterize the crack tip field. 

Close to the crack tip, the strain, e^, stress, o , effective stress, o Q , 
and effective plastic strain e P , can be expressed as 


e w = s, = t. 
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/ Ey 
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(1.29b) 
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The values of the slopes of the line segments close to the crack tip in Figs. 

1.6 to 1.9, m-j , » m^, and m^, are tabulated in Table 1.1. These values are 

very close to those given by Hutchinson [4], and Rice and Rosengren [5], except 
the plane strain, and mp, which are slightly different. The sums of the cal- 
culated values of (m-j + m^) and(m 3 + m^) vary from 0.94 to 1.04. These are close 

to unity as predicted by the HRR theory. 

Figures 1.6 - 1.9 are for n = 0.3. Similar plots for n = 0.2 are shown in 

Figs. 1.10 - 1.13. In these figures the individual points have been deleted, and 

the far field data are shown. 
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TABLE 1.1 


COMPARISON OF 

THE FEM 

AND HRR THEORY FOR THE VALUES OF m ] 

, , 

m^, find 

Stress-Strain 

State 

n 

m 

1 

m 2 


m 3 


m 4 



FEM 

HRR 

FEM 

HRR 

FEM 

HRR 

FEM 

HRR 

Plane Stress 

0.2 

0.77 

0.83 

0.167 

0.167 

0.83 

0.83 

0.167 

0.167 

Plane Stress 

0.3 

0.71 

0.77 

0.23 

0.23 

0.81 

0.77 

0.23 

0.23 

Plane Strain 

0.2 

0.78 

0.83 

0.225 

0.167 

0.83 

0.83 

0.167 

0.167 

Plane Strain 

0.3 

0.75 

0.77 

0.29 

0.23 

0.79 

0.77 

0.23 

0.23 
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Figures 1.14 to 1.16 plot the product of • ■£&-) at various 
distances from the crack tip,r, as a function of the Y produ^t of the far 
field parameters, • ^-). In Figs. 1.17 to 1.19 , the normalized 

Oy e u 

J-integral values are plotted against (a a tja^e^) and deformation work 
density Wp respectively for all the cases analyzed. It is demonstrated in 
these figures that both in the small scale yielding region (a^c^/a^c^ * 0.5) 
and in the general yielding region {ccoCoo/oy ty > 1.4) the relations between 
( 0 yy e yy/ a Y e Y ) an ^ (o°°e®/ 0 y£y)» (J/a^a) and e^), and (J/oyeya) and 

(Wp/ayEy) appear linear. There is a nonlinear region in between (e.g. 
O.5<(0 oo eo O /oyEy )<1 .4). Several correlations between J and far field parameters 
in this region have been suggested in the recent years [21-24] 

Guided by Eqns. (1.20a) and (1.20b), an obvious modification to approxi- 
mate the relation between J and (ccEc) for the single edge cracked specimen 
in tension, both plane stress and plane strain is [7] 
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E )-l]F(^ *n) ; o Q 


> o v (1.31b) 


Y Y 


where F(^,n) is the factor similar to g(^,n) and f(^,n) in Eqns. (1.12) and 
(1.13) for pure power hardening materials. We assume that F is related to g 
and f in Eqns. (1.12) and (1.13) 

F(J,n) = k(-y] 1+n f(^,n) for plane strain (1.32a), 
F (|,n) = kg(^,n) for plane stress (1.32b). 


The values of f and g are tabulated as f 5 and g-| in the References [6] and [7] 
for pure hardening materials with center cracked strip. These values are used 
for the correlation for this study. The value of k = 1.12 was used to fit 
the data. 

According to the results of Kumar et al. [25] for Ramberg-Osgood power 
hardening materials, the value of k varies with a/w and n. The value of F, 
at a given applied stress, must be the asymptotic value as a/w approaches 
zero, and J will be linearly proportional to a. 
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Figure 1.17 shows the comparison of the results from FEM calculation 
and the results based on the simple estimates of Eqn. (1.31) for n = 0.2, 
a/w = 0.1. The results of the FEM calculations and the estimates by Eqn. (1.31) 
are listed in Tables 1.2 and 1.3 for comparison. The maximum error is 
less than 92 for both the plane stress and plane strain cases at the point of 

Ocoeco/oyey =1. 

Figures 1.17, 1.18, and 1.19 show the correlation between J and 
and W D . These results substantiate that J can be linearly related to these 
two far field parameters in the general yield region, i.e. o^oy 

J = [n^e* + n^yEyla (1.33) 

J = [nj'W D + T, 2°y £ y^ a O’ 34 ) 

• I* II 

^2 and ^ are the intercepts of the lines when and W D are zero. 
is nearly equal to zero, if the strain hardening exponent n is not too small. 
Both and W D are related to deformation work done to a solid, DW. There- 
fore Eqns. (1.33 and 1.34) can be written in the form 

J = [n 1 DW + n 2 o y c y ]a (1.35) 

where 


or 



and it is the primary driving force for crack extension. 
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II. J-1NTEGRAL AND FATIGUE CRACK GROWTH 

In strain controlled fatigue or low cycle fatigue, cracks are initiated and 
the cracks propagate in cyclic general-yielding loading.. Therefore, the LEFM 
can no longer be used to analyze fatigue crack growth rate. 

Dowling and Begley [26] and Dowling [22] suggested to use AJ to correlate 
with fatigue crack growth rate. Figure 2.1. The value of J was evaluated as 
the "deformation work difference" for through cracks in plates. Dowling [22] 
has developed an equation for surface crack by combining the FEM calculations 
of center-cracked and edge-cracked panels by Shih and Hutchinson [7]. Haddad 
and Mukerjee [27] and Tanaka, Hoshide, and Nakata [28] followed the same pro- 
cedure to evaluate J and correlated J with da/dN. Kaisand and Mowbray [21] 
correlated fatigue crack growth in general yielding with AJ. AJ was divided 
into two parts: elastic and plastic. 


AJ = AJ e + AJp (2.1) 

Using Shih and Hutchinson's calculation and following a procedure similar to 
Dowling's, AJ p for a surface crack is approximated by 

AJ p = 1.96/T7n' AW“a (2.2) 


where n' = cyclic strain hardening exponent. AW p = plastic deformation work 
density at infinity. 

Tomkins [29] and Tomkins, Summer and Wareing [30] correlated crack growth 
rate with crack tip opening displacement, CTOD, and J. J consists of two parts 


J = 


2 

Tio a 


onroCpa 

1+n 


(2.3) 


where n is monotonic strain hardening exponent. 

Solomon [3] measured fatigue crack growth in specimens in general yielding. 
In this section, the crack tip field of Solomon's specimens will be calculated 
with FEM. The J values will be evaluated, and the measured crack growth rate 
will be correlated with the calculated J values. The crack growth data of a 
number of steels and aluminum alloys were examined and used to develop a general 
expression for the characteristic fatigue crack growth behavior in terms of aJ. 
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2.1 The Analysis of Solomon's Experiment 

The test section of the Solomon's specimen is shown in Fig. 2.2. The test 

section of the cylindrical specimen was reduced by two semi-circular notches. The 

2 

gross cross-section of the specimen was 1.24 x 12.6 rm . A sharp notch was used to 
initiate the fatigue crack. A specimen was rigidly attached to the loading frame of the 
test machine. Fatigue crack growth rates at six cyclic plastic strain ranges were 
measured. The applied strain range was controlled by an extensometer located 
at 7.62 mm from the edge of the specimen as shown in Fig. 2;2. The specimens 
were made of 1018 steel. The chemical composition is given in Table 2.1. The 
cyclic stress-strain relation of this steel is 

p- = 1360(^-) 0 - 26 (2.4) 

. Ac 

where and are cyclic stress amplitude and cyclic plastis strain am- 
plitude. MPa. This cyclic stress-strain relation was used in the FEM 

calculations. 

The total strain consists of two parts 

Ac = Ae g + ACp (2.5a) 

and 

Ae e = I 2 - (2.5b) 

The cyclic yield stress (170 MPa) is obtained as the intersection of two straight 
line segments in the log Aa/2 versus log Ae/2 plot, one in the elastic region and 
one in the plastic region (see Fig. 2.3). 

Solomon [3] plotted his crack growth data in terms of a pseudo stress in- 
tensity factor defined as 

A(PK) = E(Ac )/a (2.6) 

where Ac = Ae e + ACp 

Subsequently, the data was analyzed by Haigh and Skelton [31] in terms of 
strain intensity factor defined as 

AK e = (^- Ae e + Aep)*^ (2.7) 
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In order to compare with Eqn. (2.7), Eqn. (2.6) can be written in the 

form 

^4^- = (Ae e + Ae p )/rra (2.8) 

Solomon correlated his data with the total strain range Ae . In Haigh and 
Skelton's analysis, the elastic strain amplitude was used instead of the elastic 
strain range because of the consideration of crack closure. The data and the cor- 
relations proposed by Solomon [3] and Haigh and Skelton [31] are shown in Fig. 2.4. 
The scatter band of Solomon's correlation is a factor of slightly more than 3, and 
the scatter band is reduced to a factor of 2.5 by Haigh and Skelton correlation. 

Brown et al . [32] calculated the plane stress crack tip field of Solomon's 
specimen with the FEM. In their calculation, a constant stress boundary condition 
was used. The crack growth rate was correlated with the size of severe strain 
zone, r $ . 

A plane stress finite element method calculation was made for the Solomon's 
specimen in the present work. The mesh is shown in Fig. 2.5a with an enlarged 
view of the detailed mesh near the crack tip in Fig. 2.5c. The lateral curva- 
ture of the specimen was simulated by layers of varying thicknesses as illustrated 
in Fig. 2.5b. 

The solid lines in Fig. 2.5c delineate the meshes, and the dashed curves 
are the paths for J-integral evaluation. The details of the calculation procedure 
have been discussed earlier in Section I. All the elements in the layer closest 
to the crack line are of the same thickness, and the values of J-integral are 
evaluated along the paths within this layer. 

Solomon's specimens were tested in the strain controlled fatigue test fixture. 
Both ends of the specimen were firmly attached to a rigid test frame. Therefore, 
a constant displacement boundary condition simulates the test condition closely, 
and a constant displacement boundary condition was used in the FEM calculation. 

Crack tip stresses and strains for three crack lengths, a = 1 .27, 2.54, and 
5.08 mm, were calculated. The data on effective stress, effective plastic strain, 
stress o , strain e are plotted in Figs. 2.6 and 2.7. Figure 2.8 shows the 
normalized strain distribution near the crack tip for three different crack 
lengths. It is interesting to note that the normalized strain data near the 
crack tip are on the same straight line for all three crack lengths. The 
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results in Figs. 2.6, 2.7 and 2.8 suggest that J Integral characterizes the 
crack tip field for Solomon's specimens, and J can be used to analyze his crack 
growth data. 

The J-values were calculated following the dashed contours in Fig. 2.5c. 

The J-values of three specimens with three different crack lengths are plotted in 
Fig. 2.9 as a function of the plastic strain at the location of the extensometer, 
i.e. at x = 7.62 mm from the left hand side edge of the specimen. J is also 
plotted as a function of the total strain, e at the location of the extensometer. 

For a strain controlled fatigue test, the stabilized hysteresis loop is 
symmetrical as shown in Fig. 2.10. During the lower half of the cycle when the 
applied stress is negative, a crack will close and will lose its effectiveness, 
and the whole specimen will deform uniformly [31]. This deformation does not con- 
tribute to the characteristic crack tip field. Therefore, the applied AJ, which 
causes a crack to open, induces the characteristic crack tip field, and propagates 
the crack, corresponds to the positive part of the loading cycle, ABC in Fig. 2.10. 

During the positive part of the loading cycle, the stress-strain follows the 
curve ABC. This segment of the stress-strain curve consists of the elastic and 
plastic parts. The plastic part is ACp and the elastic part is (Ae-Ae p )/2. The 
AJ-value for Solomon's specimen was evaluated at the plastic strain range Ac p , 
and it corresponds to the AJ along the loading curve ABC in Fig. 2.10. 

The stress-strain relation along ABC differs from the cyclic stress-strain 
curve. However, the cyclic stress-strain curve was used to calculate the crack 
tip field for Solomon's test. From the figure, the AJ-values of the 18 combina- 
tions of three crack lengths and six cyclic plastic strain ranges can be obtained. 

The crack growth rates at these 18 combinations of crack lengths and At p 
are obtained from Solomon's data with A(PK) = EAe^t/a . The correlation of AJ 
with da/dN for Solomon's test is shown in Fig. 2.11. The scatter band in Fig. 2.11 
is much narrower than the plot in terms of either A(PK) or AK^ in Fig. 2.4. The 
data in Fig. 2.11 gives the empirical relation 

^- = 0.25 x 10 6 AJ K6 (mm/cycle) (2.9) 
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for Solomon's 1018 steel. Equation (2.9) is also plotted in Fig. 2.1 as the 
dashed line. The data of 1018 steel is very close to the data band of A533B 
steel . 


2.2 Fatigue Crack Growth in the Near Threshold Region 

An abundance of information on fatigue crack growth is available in the 
literature. The overwhelming majority of the data in the near threshold region 
is correlated with the stress intensity factor range, AK. Since K and J are 
related, the same data can be converted to correlation with aJ. 

Crack growth rate decreases with AK and AJ. da/dN starts to decrease 

rapidly at the transition point at da/dN = and AK = AK^. or AJ = AJ^.. Liu and 

Liu [33] have shown that, in the near threshold region, fatigue crack growth 

rate can be expressed as 

m - - tf - «t£n^ - ^ < 2 -' 0a > 

Yoder et al . [34] suggested that for steels, the transition point is related to 
the mean free path for dislocation movement, £, Young's modulus E, and cyclic 
yield strength Oy^, and we have 

Sr = ci ir ^ • & < 2 - 10b > 

where D, F, C, 6, m, and AK^ h are empirical constants. 

Since K and J are related, therefore these equations can be written as 


jj|. D'CjjJi-- B'f - ff‘ (2.»a) 

and 

m - c '* tP tsr - *'?' « 2 -' lb > 

th 

/ 

Equation (2.11b) is used to make the regression analysis of the fatigue crack 
growth data of a number of steels, which were analyzed earlier by both Yoder et al. 
[34] and Liu and Liu [33]. In the regression analysis, the values of C', B' , and 
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m' are allowed to vary. The values of C-j , B 1 , and of each data set are de- 
termined at the point of minimum error. These values are tabulated in Table 2.2. 
Both of the average values of B' and m' are close to one. Once the value of B' 
is accepted as equal to one, a second regression analysis was made to determine 
the values of C' and m', Cg and n^. These values of and m^, and B = 1 are 
used to plot the curves in Fig. 2.12 . With the average values of B^ and m^' , 

the data of all the steels in Fig. 2.12, give the empirical equation, which is 
shown as the dashed line in the same figure. 

&' 0 - 1 ’> (2 ' ,2a) 


With 


i = AK t /8TTOy£ = AJ^ E/8Ttayj. and (AK^/AK^.) = AJ t /AJ^ = 2, Eqn. (2.12a) can 


be written in the form 

da _ „ AJ th , AJ lX m 

dN °YC 4J th ’ 


(2.12b) 


The data in Fig. 12 give an average value of A equal to 1/40 tt. 

The data of nine additional steels are plotted in Fig. 2.13 according to 
Eqn. (2.11a). Table 2.3 shows the regression results of these nine steels. 

Figures 2.15 to 2.24 show individual regression curve for each fatigue crack 
growth curve. The original data and the references of the data sources are 
shown in Fig. 2.13. 

The data of several aluminum alloys are shown in Fig. 2.14. The data starts 
to deviate from Eqn. (2.11a) when da/dN exceeds 2 x 10 mm/cycle. Even for low 
cycle fatigue tests, da/dN does not often exceed this growth rate. Therefore, 
one may conclude that fatigue crack growth data in the near threshold region 
correlates well with Eqns. (2.11a and 2.11b). 

If the second term of Eqn. (2.12) is neglected for the FCG rate in the higher 

- 2 9 2 2 

AJ region, and if the relation l = AK^ /8irOyj. = AJ^. E/8irayj. , (AKj./AK^) 

= AJ^/AJ th = 2, and C' = 0.1 are used, Eqn. (2.12) becomes 


da 

dN 


0.01 


AJ 

°YC 


(2.13) 


Eqn. (2.13) agrees reasonably well with the unzipping FCG calculation obtained by 
Kuo and Liu [35,36] 


da _ 
dN ‘ 


0.02 


AJ 

°YC 


(2.14) 
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1018 Steel 
A533B Steel 


In conclusion, a general expression of the following form will fit well 
with each individual fatigue crack growth curve 


% = A[AJ - AJ th ] m (2.15) 

In the higher AJ region, the second term can be neglected and 

^ = AAJ m (2.16) 


TABLE 2.1 

COMPOSITION OF THE 1018 STEEL AND A533B PRESSURE VESSEL STEEL 
C Mn P S N 0 H Ni Mo Si 

0.18 0.89 0.013 0.022 0.0047 0.0019 <0.002 000 

0.25 1.15^.1.50 0.035 0.40 0 0 0 0.4^. 7 0.45M1.60 0.15^0.30 
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THE REGRESSION RESULTS OF YODER'S DATA BY USING da/dN VERSUS AJ 
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Average 1.005 1.0203 4.76 

Value 


REGRESSION ANALYSES OF FATIGUE CRACK GROWTH DATA 
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Table 2.3 (continued) 
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III. FATIGUE CRACK GROWTH IN GENERAL YIELDING AND LOW CYCLE FATIGUE LIFE 

Fatigue damage includes the processes of crack initiation and propagation. 
Owing to the defects and damage during material processing and fabrication, and 
owing to environmental attacks such as oxidation and hydrogen embrittlement, the 
crack initiation period might be shortened, and fatigue lives of such engineering 
components or structures are mainly fatigue crack propagation periods. 

Elastic-plastic fracture mechanics can be used to study the life of a low 
cycle fatigue specimen based on crack propagation. Kaisand and Mowbray [21] de- 
rived the expression relating the life of strain controlled fatigue specimen to 
the far field deformation work density, aW°°. The J-integral relation developed 
by Shih and Hutchinson [7] for the Ramberg-Osgood hardening materials was used. 

According to Kaisand and Mowbray [21] the relation between AW°° and fatigue 
life, N^, can be expressed as follows: 

AW" = AW e ~ + AW” 

= C 2 N f" 1/B + C l N f" 1/Y (3.1) 

where AW”is the applied deformation work density which consists of elastic and 
plastic parts. , C 2 » 8 and y are constants. AW g “ and AWp°° are related to 
the applied stress and strain ranges. 

J is related to deformation work density, J = n ode = nW D [6,7,9,18,37]. 

Dowling and Begley [26] and Dowling [22] correlated fatigue crack growth with AJ; 
and Chen, Zheng, and tsai [38] have correlated low cycle fatigue lives with the 
deformation work density. 

In this section, the cyclic life will be related to the applied stress and 
strain ranges, Ao^ and Ae^. The relation will be derived from the relation between 
J and the far field parameters together with the cyclic crack growth relation in 
Sections I and II. 
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A fatigue crack is often nucleated on a free surface in the shape of a semi- 
ellipse, and the nucleated crack is often a shear crack. For a part through shear 
crack, the shear mode of the crack tip field varies from mode II to mode III along 
the periphery of the crack. For a part through tensile crack, Fig. 3.1, a plane 
stress crack tip field exists at A, and a plane strain field exists at B. The 
variation of the character of the crack tip field is determined by the .crack orien- 
tation and the applied stresses. 

Along the periphery of a semi-elliptical shear crack, the value of J varies. 
Equations (1.33) and (1.34) have to be modified by a function of 0,<J>(0). <p(d) pre- 
scribes the 0-variation of J along the crack periphery. In the steady state growth of a 
small crack, the geometry of the growing crack can be assumed to remain similar. 

For geometric similar cracks, 4 >(e) remains unchanged, and only a single crack 
length parameter is needed to characterize the crack tip field. Therefore, equa- 
tions in the forms of (1.33) and (1.34) are applicable to the geometrically simi- 
lar semi-elliptical shear cracks, as well as to the geometrically similar semi- 
ellipitical tensile cracks. 

The relations between J and the far field parameters were developed for 

monotonic loading. For completely reversed cyclic loading, R = -1 and o mQ = 0, 

mean 

a crack will close under compression, losing its effectiveness. Therefore, the 
processes of developing the characteristic crack tip field, crack tip opening, 
and crack growth occur during the tensile half cycle, i.e. during the part of 
the loading curve ABC in Fig. 2.10. The effective aJ for crack growth corresponds 
to this part of the loading curve. Hence, the effective stress range, Ao , and 

A 00 

Ao A e 

effective strain range, Ae^, are and (Ae p + -jp). Therefore, a ^ and in Eqns. 
(1.24) and (1.25) should be replaced by these two quantities, correspondingly. The 
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quantity Oy and ey are used to normalize the expression for J in Figs. 1.17 
and 1.18. Hence, Oy£y can be changed into °yq e yc without major change in the 
form of the derived relation. A small change in the form of a constant factor 
can be absorbed into the constant 02 - 

In Section II, the J of a mode I through crack has been investigated. For 
geometrically similar part-through semi-elliptical cracks in infinite solids, 
with the e-variation omitted, the J integral can be written in the similar form. 
For a piecewise power hardening material, 


AJ = (r^ DW + r >2°YC e YC^ a 


(3.2) 


where and n 2 are functions of n and they are constants for a given material. 
DW is the deformation work parameter. It can be either Ao^Ae^ or AWp. 

In Section II, the data indicate that a general crack growth relation can 
be put in the form 


da 

dN 


A 

a YC Aj th 



(3.3) 


If the growth of semi-ellipitical cracks follows the same relation, and substi- 
tuting Fqn. (3.2) into Eqn. (3.3) and integrating, we obtain the following ex- 
pression. During the propagation of a small crack, DW is assumed to be constant. 


(Nf N q ) A(m-l)[ni DW+n 9 o vr e Yr ]^ a m-1 a f m-1 ^ * 

1 2 YC YC / 0__ . 1} (_L - l) 

a th a th 

where a and a,. are initial and final crack lengths. N and N f are crack 
of u i 

initiation life and total fatigue life. a th can be considered as the maxi- 
mum size of a non-propagating crack at the applied cyclic load. 
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a th = Aj th /(n l DW + n 2°YC e YC ) ^• 5) 


At a given value of DW, a^ is constant. 

If a^r is much larger than a Q , the second term within the bracket of (3.4) 
is much smaller and can be neglected. After rearrangement, Eqn. (3.4) becomes 


[n T DW 


+ n 2 a 


AJ 


YC e YC^ a p) 


m-1 


th 


[n, DW 


+ n 2 a Y C e Y C 


](N, 


- N o> 


YC 

A(m-1 ) 


(3.6) 


Fatigue life is a function of the deformation work, DW. Thus we have 

(N f - N 0 ) = f ( DW ) = f^Ao/eJ = f^) (3.7) * 


a is related to the microstructural feature of the crack nucleation process, 
o 

For a given material, a Q is a constant. n 2 * s negligible if the value of n is not 
too large. N Q is negligible if there exist material and fabrication defects or if 
environmental attact is severe. If a Q is much larger than a th> and if n 2 and N Q 
are negligible, then we have 

o vr AJ.. m-1 

(DW) m N f = (— = C (3.8) 

* A(m-l)n" a o 

The right hand side of Eqn. (3.8) is a constant. Or 

( Ao a, Ae J m N f = c l 

and 

A W D m N f = C 2 


(3.9) 

(3.10) 


The total life consists of fatigue crack initiation and propagation, N Q and 
(Nf - N q ). The crack initiation period is strongly affected by the size of the 
material and fabrication defects, a Q . A large a Q will reduce the total fatigue 
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life. The quantity, a^ can be considered as the size of a non-propagating 
crack, at the applied stress and strain ranges. If a Q is less than a th , the 
crack will not grow and no fatigue failure will occur. Fatigue life is also 
affected strongly by the crack growth characteristics A, m and aJ^. 

Because of the uncertain nature of the cyclic crack tip field due to the 
lack of the precise cyclic constitutive relation, the lack of the calculated 
cyclic crack tip field, the lack of the information on shear cracks and the 
uncertain nature of the frictional shear stress between mating crack surfaces, 
the relations given by Eqns. (3.4, 3.6, 3.7, 3.8, 3.9 and 3.10) may not be accu 
rate quantitatively. However, these equations give the basic forms of the ex- 
pressions for fatigue life, and how fatigue life is related to all of these 
factors. And these equations also help us to develop meaningful correlations 
of the experimental data and to assess the contributions from various physical 
and mechanical factors. 

Chen, Zheng, and Tsai [38] conducted strain controlled fatigue tests on 
copper-aluminum wrought alloys. The chemical compositions of the aluminum 
alloys are shown in Table 3.1. Their mechanical properties are tabulated in 

Table 3.2. The specimen geometry is shown in Fig. 3.2. 

Aa^ Ae 

The plots of (Aa^Ae^) vs and ^ — (Aep + ~ 2 ~) vs N^. are shown in Figs. 

3.3 and 3.4. Figure 3.4a shows the fatigue lives of both smooth speci- 

mens and precracked specimens. The specimens were precracked by electrostatic 
discharge machining. Both lines in Fig. 3.4a are parallel. The difference 
in cyclic life can be interpreted as the number of cycles to grow a crack from 
its "natural" crack nucleus size to the size of the precrack. If we denote a o -j 

and N^i as the precrack size and the fatigue life of a precracked specimen, 

Eqns. (3.8 and 3.10) give the ratio N f /N fl = constant. Therefore, the plot of 
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the fatigue lives of smooth and precracked specimens should be parallel as shown 
in Fig. 3.4a. But theiinesof the precracked specimens are shifted 

to the left. 

The plots of AWp vs are shown in Figs. 3.5 and 3.6. aW^ can be defined 
as 




■Ae p + Ae g /2 

ode 

0 


(3.11a) 


AWp is for the part of the loading curve that causes a crack to open and to 
grow, i.e. ABC in Fig. 2.10. For these specimens, the values of deformation 
work density were directly measured from stabilized cyclic stress-strain 
hysteresis loops at approximately half of the fatigue life. The plots of the 
area under the hysteresis curve vs N f are shown in Fig. 3.5. 

AWg can also be defined as 


AW 0 = 


■Ae e /2 

ode 

0 


% 

ode 


(3;'ilb) 


Take the approximation of the cyclic stress and strain ranges in the form 
p Z J YC 


n 1 Ao 

(Ao - Oy C ) = Ae p for y ~ > o yr and integrating Eqn. (3.11b), we obtain 


4W d = 8E~ + c TSf’ 4e p + 0 yc ]4t , 


sT + [ )TF' (40 - C ’vc ) + °vc :l4e p 


(3.12) 


32 



The first term for the plastic strain is much smaller and can be neglected. 

2 


aW D = 8E + °YC A£ p 


(3.13) 


Substituting Eqn. (3.13) into Eqn. (3.8), we obtain 

,m-l 


C « + °YC ie p^ N f = 


m 


8E 


A(m 


°YC Mth\ 
•D(nj) m \ a o / 


(3.14) 


where a Q is the initial crack length. The plots of AWp = [g|- + o^Ae^] vs N^. 
are shown in Fig. 3.6. 

Straight lines correlate the data very well in all of these figures. Figs. 

5 

3.3 to 3.6, over the fatigue life range from 10 to 10 cycles. Therefore Eqns. 
(3.9) and (3.10) are valid over three to four orders of magnitude of fatigue 
life. Additionals data are necessary to distinguish these four correlations. 
Equations (3.9) and (3.10) can be written in the forms 


E(- 


Ae 


1 


and 


e\2 w m . Ao , M i 

2~) N f + r Ae p N f 


■ i i 

E/_j?\2 w m . N m 
22 N f + °YC Ae p N f 


= C' 


= C" 


(3.15a) 


(3.15b) 


The left handsidesof these two equations consist of the elastic and the plastic 
parts. Each of these parts is similar to the basic form of the elastic and 
plastic contributions of strain controlled fatigue, i.e. 


and 


Ac e N f - C 1 


y 2 

Ac p N f = C 2 


(3.16) 


When the elastic part is negligible, the plastic part becomes the original 
Coffin-Manson low cycle fatigue law. 
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Table 3.1 Chemical Composition for 
the Alloys Tested 



Cu 

Mn 

Si 

Fe 

Mg 

Zn 

Ni 

A1 

LD-10 

4.56 

0.86 

1.01 

KB 

0.70 


0 

92.631% 

LY-12 

3. 8*3. 9 

0.3*0. 9 

0.5 

0.50 

1.2*1. 8 

0.3 

0.1 

91*93.31% 


Table 3.2 The Material Properties for 
LD-10 and LY-12 Aluminum Alloys 




Oy (MPa) 

° u (MPa) 

E(MPa) 

o Y( . (MPa) 

Al. Alloy 

Heat treatment 
Condi tion 

Yield 

Strength 

Ultimate 

Strength 

Young's 
Modu 1 us 

Cyclic Yield 
Strength 

LD-10 

Cold rolling 

281 

439 


323 

LY-12 

Natural age 

416 

578 

7. 35x1 O 4 

465 

■ mm 

Artificial age 

313 

411 

7.45xl0 4 

274 
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CONCLUSIONS 


1. A characteristic crack tip field exists for a small through crack in a 
large plate in both plane strain and plane stress conditions. The crack 
tip field can be characterized by J. 

2. For small cracks, J is linearly proportional to crack length and propor- 
tional to the far field parameters: the product and the applied 

deformation work density. 

3. Fatigue crack growth rate in general yielding correlates well with AJ. 

4. An expression for the general fatigue crack growth behavior is developed 
in terms of aJ. 

5. An equation for fatigue life is derived from the general crack growth be 
havior. This equation relates fatigue life to various physical and mech 
anical aspects of fatigue damage, such as crack nuclei size, fracture 
toughness and fatigue crack growth behavior. 

6. Low cycle fatigue lives of two aluminum alloys correlate well with de- 
formation work density. 

7. The derived equation for fatigue life is reduced to Coffin-Manson low 
cycle fatigue law in the high strain region. 


35 



ACKNOWLEDGEMENT 

The work was conducted at the George Sachs Fracture and Fatigue Research 
Laboratory at Syracuse University. The financial support by NASA, Grant No. 
NAG 3-348 and the Visiting Scholarship by the People's Republic of China to 
Mr. Zheng Minzhong are acknowledged. The Computer Center of the Syracuse 
University for making computer time available to the authors and the Hibbit, 
Karlsson, and Sorensen, Inc. for the privilege of using the ABAQUS finite 
element program are also acknowledged. The assistance by Mrs. Helen Turner, 

Mr. Rolf Zeimer and Mr. Qi Chen in manuscript preparation is very much 

appreciated. 


36 


REFERENCES 


[1] L. C. Coffin, Jr., "A Study of the effects of cyclic thermal stresses 
on a ductile metal," Trans. ASME 76, 1954, p. 931. 

[2] R. W. Smith, M. H. Hirschberg, and S. S. Manson, "Fatigue behavior of 
materials under strain cycling in low and intermediate life range," 

NASA TN D-1574, Nat. Aero, and Space Adm., 1963. 

[3] H. D. Solomon, "Low cycle fatigue crack propagation in 1018 steel," 

Journal of Materials, JMLSA, Vol. No. 3, September 1972, pp. 299-306. 

[4] J. W. Hutchinson, "Singular behavior at the end of a tensile crack in a 
hardening material," J. of Mech. and Phys. of Solids, Vol. lj>, 1968, p. 13. 

[5] J. R. Rice and G. F. Rosengren, "Plane strain deformation near a crack 
top in a power-law hardening material," J. of Mech. and Phys. of Solids, 
Vol. 16, 1968, p. 1. 

[6] N. L. Goldman and J. W. Hutchinson, "Fully plastic crack problems: The 
center-cracked strip under plane strain," Int. J. Solids Structures 11_, 
1975, p. 575. 

[7] C. F. Shi h and J. W. Hutchinson, "Fully plastic solutions and large scale 
yielding estimations for plane stress crack problems," J. Eng. Mat. and 
Tech., 1976, pp. 289-295. 

[8] J. C. Amazigo, "Fully plastic crack in an infinite body under anti-plane 
shear," Int. J. Solids Structure 10, 1974, p. 1003. 

[9] J. R. Rice, Journal of Applied Mechanics, Transactions, American Society 
of Mechanical Engineers, Vol. 35, 1968, pp. 379-386. 

[10] C. F. Shih, "J-integral estimates for strain hardening materials in anti- 
plane shear using fully plastic solution," Mechanics of Crack Growth, 

ASTH STP 590, 1976, pp. 3-26. 


[11] A. A. Ilyushin, "The theory of small elastic-plastic deformation," 
Prikadnaia Matematika Mekhanika, P.M.M. 10, 1946, p. 347. 

[12] L. M. Kachanov, "Fundamentals of the theory of plasticity," translated 
from the Russian by M. Konyaeva, Mir Publishers, Moscow, 1974. 

[13] J. W. Hutchinson, "Singular behaviour at the end of a tensile crack in a 
hardening material, and plastic stress and strain fields at a crack tip," 
J. Mech. Phy. Solids 1J5, 1_3, 1968, 337. 

[14] R. M. McMeeking and D. M. Parks, "On criteria for J-dominance of crack 
tip fields in large scale yielding," Elastic-Plastic Fracture, ASTM, STP 
668 , Am. Soc. for Testing and Materials, 1979, pp. 175-194. 


37 



[15] C. F. Shih and M. D. German, "Requirements for a one parameter charac- 
terization of crack tip fields by the HRR Singularity," Int. J. of 
Fracture 1_7, No. 1, Feb. 1981, pp. 27-43. 

[16] W. L. Hu, "A finite element study on crack tip deformation," Ph.D. 
Dissertation, Solid State Science and Technology, Syracuse University, 1976. 

[17] W. L. Hu and H. W. Liu, "Characterization of crack tip stresses and 
strains — small scale yielding and general yielding," Proceedings of the 
Second Inti. Conf. on Mechanical Behavior, Aug. 16-20, 1976, Boston, Mass. 

[18] H. Tada, P. Paris, and G. Irwin, "The Stress Analysis of Cracks Handbook," 
published by the Del Research Corporation, 1973. 

[19] A. S. Kobayashi , "Numerical analysis in fracture mechanics," The Office 
of Naval Research under Contract N00014-76-C-0060 NR 064-478. Technical 
Report No. UWA/DME/TR-83/45. 

[20] Zhaung Tao, Ph.D. Dissertation, Dept, of Chemical Eng. and Mat. Sci., 
Syracuse University, 1984. To be published. 

[21] L. R. Kaisand and D. F. Mowbray, "Relationships between low-cycle fatigue 
and fatigue crack growth rate properties," Journal of Testing and Evalua- 
tion, JTEVA, Vol . 7, No. 5, Sept. 1979, pp. 270-280. 

[22] N. W. Dowling, "Crack growth during low-cycle fatigue of smooth axial 
specimens," Cyclic Stress-Strain and Plastic Deformation Aspects of Fatigue 
Crack Growth, ASTM STP 637, 1977, pp . 97-121. 

[23] R. P. Skelton, "Growth of short cracks during high strain fatigue and 
thermal cycling," Low-Cycle Fatigue and Life Prediction, ASTM STP 770 , 1982, 
pp. 337-381 . 

[24] V. M. Radhakrishnam, "Damage accumulation and fracture life in high- 
temperature low-cycle fatigue," Low-Cycle Fatigue and Life Prediction, 

ASTM STP 770, 1982, pp. 135-151. 


[25] V. Kumar, M. D. German, and C. F. Shih, "An engineering approach for 
elastic-plastic fracture analysis" Electric Power Research Institute, 
NP-1931 Research Project 1237-1, Topical Report, July 1981. 

[26] N. E. Dowling and J. A. Begley, "Fatigue crack growth during gross plas- 
ticity and the J-integral ," Mechanics of Crack Growth, ASTM STP 590 , 1976, 
pp. 80-103. 

[27] M. H. El Haddad and B. Mukherjee, "Elastic-plastic fracture mechanics 
analysis of fatigue crack growth," Elastic-Plastic Fracture: Second Sym- 
posium, Vol. II - Fracture Resistance Curves and Engineering Applications, 
ASTM STP 803, 1983, pp. 1 1 -689- 1 1 - 707 . 


[28] K. Tanaka, T. Hoshide, and M. Nakata, "Elastic-plastic crack propagation 
under high cyclic stresses," Elastic-Plastic Fracture: Second Symposium, 
Vol. II - Fracture Resistance Curves and Engineering Applications, ASTM 
STP 803, 1983, pp. 11-689-707. 


38 



[29] B. Tomkins, J. of Engineering Materials and Technology, Transactions of 
ASME, Series H. , Vol. 97, 1975, pp. 239-297. 

[30] B. Tomkins, 6. Sumner, and J. Wareing, in Inter. Symposium on Low Cycle 
Fatigue Strength and Elasto-Plastic Behaviour of Materials, Deutscher 
Verband fur Material Prufung, Berlin, 1979, pp. 495-507. 

[31] J. R. Haigh and R. P. Skelton, Materials Science and Engineering 36, 

1978, pp. 17-25. 

[32] M. W. Brown, H. W. Liu, A. P. Kfouri , and K. J. Miller, "An analysis of 
fatigue crack growth under yielding conditions," Advances in Fracture 
Research edited by D. Francois et al., Pergamon Press, Oxford and New 
York, 1980. 

[33] H. W. Liu and Dai Liu, "A quantitative analysis of structure sensitive 
fatigue crack growth in steels," Scripta Metallurgica 1[8, 1984, pp. 7-12. 

[34] G. R. Yoder, L. A. Cooley, and T. W. Crooker, Naval Research Laboratory 
Memorandum. Report 4576, July 1981. 

[35] A. S. Kuo, and H- W. Liu, "An Analysis of Unzipping Model for Fatigue 
Crack Growth," Scripta Met. 10. » P- 723, 1976. 

[36] C. Y. Yang and H. W. Liu, "The application of the unzipping model of 
fatigue crack growth to micro-cracks," Scripta Metallurgica 14_, July 1980, 
pp. 785-790. 

[37] C. K. Tsai, "J-integral analysis," J. of Acta Metallurgica ±2, 1(1976), p. 45. 

[33] X. X. Chen, M. Z. Zheng, and C. K. Tsai, "Study of strain fatigue in LD-10 
alloy," J. of Acta Mechanica Sinica, China 5, 1981. 


39 



J=/ p Wdy- T --f£ ds 

w =rv<-,i 

(a) Line integral contour. 

Fig. (1.1) J-Integral of Rice. 


Clamped: W= constant 



Fig. (1.2) (a) 


Infinitely wide slab 


(b) Stress-strain 
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(a) Initial mesh 
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(b) Mesh after crack blunting 
Fig. (1.5) Crack tip blunting. 
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ig. (1.7) Normalized 
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Fig. (1.8) 
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Fig. (1.10) Normalized effecti 






Fig. (1.11) Normalized stress and strain distribution ahead of the crack tip. 
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Fig. (1.15) Near field parameter a yy £ yy/ a Y e Y versus far field parameter 











Fig. (1.18) Correlation of dimensionless J versus o^e^/oyCy for single edge cracked panel 
in tension with n = 0.3 and a/w = 0.1. 
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Fig. (1.19) Correlation of dimensionless J. versus Wp/oyey for single edge cracked panel 
in tension. 
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Extensometer 
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Fig. (2.2) Dimensions of Solomon's specimen. 
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Fig. (2.4) Comparison of Egns. (2.7) and (2.8) using Solomon's data. 
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Fig. (2.8) 
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Fig. (2.9) J-integral dependence on applied strain. 










Fig. (2.10) Cyclic s tress -starin loop and the 
deformation work density. 
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Fig. (2.11) Correlation of da/dN versus AJ for Solomon's tests 
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Fig. (2.12) Fatigue crack growth data for steels, 
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Fig. (2.13) Fatigue crack growth data for steels. 
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Fig. (2.14) Fatigue crack growth data for aluminum alloys. 
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Fig. (2.15) Fatigue crack growth. data. 
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Fig. (2.17) Fatigue crack growth data. 
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Fig. (2.19) Fatigue crack growth data. 
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Fig. (2.22) Fatigue crack growth data. 


8 


f 








X3 ho 



Fig. (2.24) Fatigue crack growth data. 
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Fig. (3.3) Semi -elliptical surface crack. 
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(3.3a) Correlation of [ha 
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Fig. (3.3b) Correlation of (Ao^J versus (2N f ) for U -12 aluminum alloy. 








LD-10 aluminum alloy. 









Fig. (3.6a) Correlation of (pf- + a Yr Ae ) versus (2N f ) for LD-10 aluminum alloy 
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deformation work density as depicted by the derived theory. The general relation 
is reduced to Coffin-Manson low cycle fatigue law in the high strain region. 
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